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The Liiders postulate is reviewed and implications for quantum algorithms are discussed. A search 
algorithm for an unstructured database is described. 
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I. INTRODUCTION 

In this paper the Liiders postulate^] is used to de- 
velop a new quantum computer algorithm to search an 
unstructured database with exactly one marked record. 
The Liiders postulate was introduced as a modification 
of the original measurement theory of quantum mechan- 
ics as presented by von Neumann |2|. It describes unam- 
biguously the measurement process of observables with 
a degenerate spectrum. Our approach to the unstruc- 
tured database search differs from Grover'sQ original 
algorithm due to our focus on the details of the mea- 
surement process. This paper emphasizes conceptual is- 
sues, while implementability will be addressed in detail 
elsewhere. 

The paper is organised as follows. In Section ITU some 
notational issues and other preliminaries are covered. In 
Section ILTTI the Liiders postulate is reviewed and with the 
help of a simple example some implications of the postu- 
late are clarified. In Section ITV1 a unstructured database 
search algorithm is described. In the last section some 
concluding remarks are added. 



II. PRELIMINARIES 

Let us begin by presenting the relevant notation. We 
define, following the conventions of the field, some terms. 
A qubit represents the superposition of a pair of orthog- 
onal quantum states. The pair of states will be denoted 
as |0) and |1). Implementations of qubits have been dis- 
cussed in various settings, such as spin-i particles, po- 
larised photons, et cetera. For reasons of simplicity we 
assume that the qubits are either implemented as spin-^ 
particles or as polarised photons. Other potential imple- 
mentations of qubits are equally valid and do not change 
the basic nature of the algorithm. Spin-up and spin- 
down particles, or horizontally and vertically polarised 
photons, eigenstates with respect to a particular direc- 
tion, are written respectively as |0) and |1). 

Extensions of single qubits are obtained as usual by 
constructing a multi-particle tensor product. The ba- 
sis of the n qubit states can be written in the form: 
|e n ) <8> |e„_i) <g> • • • <g> |ei), where each |e<) is either |0) or 
For reasons of brevity we often write |e„e n _i • • • e\) 
for the n qubit. We can also write qubits in vector nota- 
tion, where each of the A states |000 • • • 00) to | 111 • • ■ 11) 
corresponds to a vector e*j of length A containing one 



nonzero component of value 1 in the j-th place. 

The following additional notation will also be used. 
The records of a database with A~ = 2 n+1 elements are 
encoded in the set D = {\oJi)} (i = 0,1, 2,..., A — 1), 
where is given by the binary representation of 

length n + 1 of the numbers to A — 1. Therefore, 
each state \u>i) is represented by a particular sequence 
|e„e n _i • • • ei), where each is either or 1. The marked 
record, in particular, will be denoted \u>k)- Normalisation 
factors associated with combination of states are con- 
sistently ignored until the end of Section IV. 



III. THE LUDERS POSTULATE 

The Liiders postulate [| & Q| describes the measure- 
ment process of observables with a degenerate spectrum, 
and it has become part of the standard canon of quan- 
tum mechanics. In the case of operators with a degen- 
erate spectrum it postulates that the projection of the 
initial wave function is onto exactly one point in each 
degenerate subspace. The point chosen is the element 
of the degenerate subspace 'closest' - in terms of tran- 
sition probability - to the initial wave function. This 
'refinement' of von Neumann's projection postulates^ 
seems reasonable, since Liiders' postulate produces mea- 
surements that disturb the wave function minimally. 

The mathematical formulation of the postulate is given 
next using standard Dirac notation. We define the nor- 
malized eigenfunctions of the observable O with K dif- 
ferent eigenvalues, each having the degeneracy dk, to be 



(1) 



where k 



1,2,. 



,K and j = 1, 2, dk- The eigen- 
functions allow the definition of the following set of K 
projection operators^ 



Pk=Y^ \*Pk,j)(*Pk,j\- 

3=1 



(2) 



A measurement of an arbitrary pure state \(j>) 5] now 
gives according to the Liiders postulate the 'reduction' 
to the following states 



Prob[0 = X k ]- 1/2 P k \ 



(3) 



with the probabilities for the distinct eigenvalues of 

Prob[<3 = A fe ] - {4>\Pk\4>)- (4) 



2 



The impact of the Liiders postulate on the distin- 
guishability of similar observables is next presented. The 
aim is to distinguish two known observables with the help 
of a measurement of an input wave function. We con- 
sider observables that measure individual spin-i parti- 
cles. The observable O is either chosen to be the identity 
operator 

/-(£!), (5) 

or the operator 

that associates the eigenvalue 1 to the cigenstate spin-up 
and 1 + S to the eigenstate spin-down. 

A measurement of either the observable / or J for an 
input wave function in equal superposition of spin-up and 
spin-down, i.e. 1/V2(|1) + |0)), is carried out next. It 
gives for the first observable I a direct projection of the 
wave function onto itself. For the second observable J 
the measurement outcome is a mixed state with equal 
probability in the state spin-up and spin-down as long as 
6 is nonzero. 

The unique outcome l/y/2(\l) + |0)) for the first ob- 
servable can be distinguished simply from the mixed state 
outcome by standard interference techniques. One can 
for example, in an additional apparatus, measure the 
probability of the wave function in an appropriate basis 
like l/\/2(|l) -I- 10)) and l/%/2(|l) - |0)). In the first case 
the outcome will always be 1 / a/2(| 1) + |0)). For the sec- 
ond case, the mixed state, the probability for each of the 
basis states is 1/2. Therefore, the ability to distinguish 
the two observables below any chosen error threshold e 
is possible, if sufficient identical copies of the system are 
prepared. Each copy available decreases the probability 
of a mistake by a factor of 1/2. If m copies are prepared, 
the probability of an incorrect choice is 2~ m . 

In effect, this simple example already captures the 
essence of the paper. Namely, an infinitesimal deforma- 
tion of an observable changes a degenerate into nonde- 
generate spectrum and can lead, as in the case described, 
to an observable difference. 

In the coming section we link the form of the observ- 
able through the use of an oracle to the location of the 
marked state. It will turn out that these different observ- 
ables, i.e. different locations of the marked state, can be 
distinguished efficiently with the help of measurements 
of specially prepared wave functions. As an aside, the 
space of Hermitian operators possesses a natural Finsle- 
rian metrical permitting a more comprehensive study of 
their properties. 

IV. SKETCH OF THE ALGORITHM 

After having introduced the necessary notation and the 
underlying principle behind the algorithm in the previous 



two sections, we can describe the process for finding the 
marked state in an unstructured database. 

Let us begin by presenting an outline of the algorithm. 
The algorithm divides the database search into smaller 
pieces to make it more tractable. We do so, initially, by 
dividing the whole set of records into two equal subsets. 
The algorithm determines within certain error bounds, 
which one of the two subsets contains the marked state. 
This is done by linking the form of the observable to 
the location of the marked state. The different observ- 
ables, i.e. different locations of the marked state, are 
used to 'measure' (following Liiders' postulate) specially 
prepared wave functions. This measurement, i.e. 'col- 
lapse' into the relevant eigenstates, gives us information 
to determine (with the help of a simple further measure- 
ment) up to a certain accuracy the presence or absence 
of the marked state in the subset under consideration. 
Repetitions of the process decrease the error probability. 
Once the set containing the marked state has been iden- 
tified to a sufficient level of accuracy, we further divide 
the set into two subsets of equal size and restart the pro- 
cess. Without loss of generality we may set N = 2 n+1 . 
The division of the records into subsets has to be carried 
out n+1 times, until we are finally left with the unique 
marked record. We define a cycle to be the process of 
halving the number of states under consideration. This 
section describes the first cycle of the algorithm. 

We present in the following three elements of the algo- 
rithm. The wave function to be measured is constructed 
first. A description of the possible observables comes 
second. Finally, we let the observable act on the wave 
function. 

A. Wave function. We begin by transforming the n 
qubit wave function of the form 1 000 • • • 00) into the 'in- 
put' wave function that will be measured by the observ- 
able. This is done by creating the superposition of all 2 n 
states of the form |000 • • • 00) + 1000 • • • 01) + |000 • • • 10) + 
■■• + |111---11) from the starting wave function by a se- 
quence of Walsh-Hadamard transforms on the n individ- 
ual qubitsQ. As noted above, we shall ignore for sim- 
plicity the normalization factors. This will not affect the 
argument. 

B. Observable. Next, we construct an observable that 
depends on the position of the marked record. It will turn 
out to be the symmetrised product of two matrices. One 
matrix is a fixed Hermitian matrix containing two block 
diagonal submatrices with the associated real, nonzero 
and unequal eigenvalues a\ and 02 (e.g. a\ = 1 + S and 
ai = 1), where the first subspace is 2-dimensional. Two 
eigenvectors, which span the 2-dimensional subspace of 
the first submatrix, are chosen to be the sum of all the 
individual states, |000 • • • 00) + 1000 • • • 01) + 1000 • • • 10) + 
■ • • + I 111 • • • 11), and the sum of all the states with al- 
ternating sign, |000 • • • 00) - |000 • • • 01) + |000 • • • 10) - 

|000 • - -11) V I 111 • • • 10) - I 111 • • -ll).In matrix form 

this can be written, due to the Schur decomposition, as 
the product of three matrices of the form A = WGR 
with G equal to 
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and R a simple unitary basis rotation matrix transform- 
ing the initial basis e\, e%, ejv into the new basis with 
the first two basis elements |000 • • • 00) + |000 • • • 01) + 
|000 • • - 10) + ••• + |111 - • - 11) and |000 • • • 00) - 
jooo • ■ ■ 01) + |000 ■ • ■ 10) - |000 • • ■ 11) ■ • • + |111 • ■ • 10) - 
I 111 - - - 11) , which are identical to the two basis elements 
described above. The rest of the new basis elements that 
make up the rows of R can be chosen arbitrarily as long 
as the resulting matrix is unitary. The product of the 
three matrices WGR is Hermitian. 

The other matrix B associates a phase of e l7r , if the 
state is the marked state and otherwise leaves the state 
unchanged. This operation is the standard delta-function 
oracle of the form fk{oJi) = Sik, implemented in the quan- 
tum mechanical context as 

\ui) - (-l) 5 **k), (8) 

where \ujk) is the marked state. Note that this transfor- 
mation is identical to the conventional oracle employed in 
the unstructured database search algorithm[?J. Feasibil- 
ity and other issues related to oracles have been discussed 
in the literature - see Nielsen and Chuang 6] for details 
and references. 

Each of the two operators A and B is Hermitian on its 
own, but the product might not, because the matrices do 
not necessarily commute. Instead we create a Hermitian 
observable by symmetrizing the product 

6= (AB + B-A)/2. (9) 

The observable C will be used in the measurement pro- 
cess. 

C. Measurement. We carry out the necessary mea- 
surements to distinguish between the different possible 
observables with sufficient accuracy. We let the observ- 
able C act on the 'input' wave function |000 • • • 00) + 
|000 • - ■ 01) + |000 • - ■ 10) + ••• + |111 - • - 11). There are 
two cases to be studied. In the first case there is no 
solution in the subset considered. The 'input' wave func- 
tion, i.e. an eigenstate of the first subspace of A, is then 
also the 'output' wave function. The operators A and C 
are identical and the result of the measurement is a pure 
state. 

The case where there is a solution in the subset con- 
sidered is more interesting to analyze. Let us define the 
vectors u\ and U2 to be a basis of the two-dimensional de- 
generate subspace of A. They are chosen to be of the form 
|000 • ■ ■ 00) + |000 • ■ • 10) + • • • + |111 ■ • -00) + Jill • • • 10), 



i.e. the sum of only the 'even' states, and 1 000 ■ • -01) + 
|000 • • - 11) + • •• + |111 • • • 01) + |111 • • • 11), i.e. the sum 
of only the 'odd' states. For the marked state in any odd 
position we have Cu\ = a\U\ and Cv,2 — a\U2 + V2 with 
V2 a nonzero vector and element of the larger degenerate 
subspace of A. For the marked state in any even position 
we have C«2 = «i'"*2 and Cu\ = axui + v% with v\ similar 
to #2 HI- Therefore, only one of the two eigenstates of the 
2-dimensional subspace of B will remain an eigenstate of 
C . As a consequence the output has to be a mixed state, 
since the 'input' wave function is equal to u\ + 

The same result can be derived by decomposing the 
operator B into the identity matrix and one isolated el- 
ement on the diagonal. The isolated element on the di- 
agonal times A will produce a rank-1 matrix. A rank-1 
matrix is also produced, if one switches the order of the 
matrices. The kernel of each of these two rank-1 matri- 
ces is (2" — l)-dimensional. Therefore, except for a 2- 
dimensional space, the eigenvalue and eigenvector struc- 
ture of the rest of A are unaffected by the multiplication 
with B and the symmetrization. The operators C and 
A possess the same eigenvalue and eigenvector structure 
except for this 2-dimensional subspace, which cuts into 
both of the degenerate subspaces of A. 

The old basis is next transformed into a new ba- 
sis to more easily distinguish the pure state from the 
mixed state. We choose our new basis in such a way 

that the wave function |000 • • • 00) + |000 •••01) H h 

1 1 1 1 - - - 10) + 1 1 1 1 - • - 11) is mapped into the new basis 
element |000 • • - 00) . The rest of the old basis can be 
mapped into any new basis, as long as the transforma- 
tion is unitary. 

Next, one measures the output state qubit by qubit in 
the new basis. There are two possibilities. In the pure 
state case all the measurements without fail produce |0), 
since one of the eigenstates of C is exactly the 'input' 
wave function in the new basis 1 000 • • • 00). 

As mentioned, in the mixed state case only one 
of the eigenstates |000 • • • 00) + |000 • • • 10) + ••• + 
|111 • - - 00) + |111 • • • 10) and |000 • • • 01) + |000 • • • 11) + 

h |111 • • -01) + I 111 - - - 11) that form a basis of the 

first degenerate subspace of A remains an eigenstate of 
C. The transition probability of the 'input' wave func- 
tion is 1/2 to both of these basis elements. Therefore, 
the separate measurements of all the qubits in the new 
basis results in at least one of the n qubits to be |1) with 
probability of not less than 1/2. Let us explain this in 
more detail. We know that one of the eigenstates has a 
transition probability of exactly 1/2 to the 'input' wave 
function. As a consequence the rest of the eigenstates 
have the same transition probability as a sum. The prob- 
ability must be at most 1/2 that the new basis element 
1 000 • • • 00) will be the outcome of the measurements and 
all the qubits are measured to be |0). If we repeat the 
cycle m times, the error probability is bounded above by 
2~ m . 

The measurement of each qubit is carried out with the 
help of a set of n properly aligned Stern-Gerlach appa- 
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ratus. Each qubit of the whole wave function will be 
measured separately. Such a splitting is achieved with- 
out damaging the coherence of the wave function (see, 
e.g., Feynman et al. for an introductory discussion of 
Stern-Gerlach experiments). 

Alternatively, if horizontally and vertically polarized 
photons are used, then the separation can be made by a 
polarising beam splitter with different reflection proba- 
bilities for horizontally and vertically polarisaed photons; 
for example, transmitting horizontally polarised and re- 
flecting vertically polarised photons. 

In the remaining paragraphs of this section we con- 
clude the description of the algorithm. After completing 
one cycle, one moves on to the next cycle. In total log 2 N 
cycles have to be performed. As one moves from cycle 
to cycle the number of qubits needed to enumerate the 
remaining states decreases one by one. Once the last cy- 
cle is finished, one has established within certain error 
bounds where the marked record can be found. 

Next, we will demonstrate that the accumulated error 
probability, inherent in moving from cycle to cycle and 
choosing smaller and smaller subsets, can be made suffi- 
ciently small. To show this, define the cumulative proba- 
bility of not choosing an incorrect subset in any one of the 
cycles to be 1-e = (l-ejv)(l-cjv/2)(l — ejv/4) 1 • 1 (!- £ 2), 
where each of the corresponds to the error probability 
in one cycle of the algorithm starting out with a fixed 
number of states, i.e., N, N/2,- ■ ■ . It is easiest, if one 
chooses a large enough number of iterations, e.g. m,for 
each cycle so that the individual cycle error, decreasing 
exponentially with m, i.e. 2" m + , is sufficiently small. 
This allows us to state a sequence of inequalities connect- 
ing the cumulative error probability with the sum of the 
error probabilities for each cycle, and then with the error 
probability for the first cycle, which in itself is bounded 
above, i.e. e < €n + £n/2 + • • • + £2 < fog2(^0 e JV < 
log 2 (Af)2 -m+1 . This bound is sufficiently small for rea- 
sonable mQ- 

Normalization factors can be introduced in the ap- 
propriate places, but would change nothing in the basic 
structure of the algorithm, since relative probabilities be- 
tween different components of the wave function are left 
unchanged. 

V. CONCLUSION 

In this section some concluding remarks are added to 
round of the paper. A number of issues can be raised in 
connection with the earlier sections: 

o Is the Liiders postulate correct [Tpj . i.e. to what ex- 
tent has it been verified experimentally? 

o What observables can one construct [rj|? How does 
one construct them? 



o Are the standard rules of non-relativistic quantum 
mechanics applicable up to the accuracy required for the 
implementation of the algorithm? There have been nu- 
merous speculations on the limits of the rules of quan- 
tum mechanics. Maybe this paper, together with the 
varied developments of quantum computing and infor- 
mation theory, will give added impetus to study these 
issues experimentally. 

o In many versions of the stochastic quantum mechan- 
ics approach there is a natural 'collapse' of any input 
wave function into the energy eigenstates of the Hamilto- 
nian. Identifying the Hamiltonian with the observable C 
is arguably the most promising approach for implement- 
ing the algorithm. What evidence is there for stochastic 
quantum mechanics? 

The case of having more than one marked state, in 
itself an interesting problem, is here only briefly com- 
mented on. Naturally, the algorithm can be modified 
to handle an unstructured database with a number of 
marked states [l^. This extension has direct applications 
to NP-complete problems. 

The purpose of this paper was to show that the Liiders 
postulate has interesting consequences for quantum algo- 
rithms, and to sketch an implementation. Issues related 
to implementability and computational complexity, i.e. 
what does a polynomial number of steps mean, will be 
discussed in a more detailed subsequent paper. 

The author wishes to express his gratitude to Dorje 
Brody for stimulating discussions. 
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[1] One can in the same way encode the possible solutions [2] In general, but not with absolute consistency, we call a 
of a NP-complete problem like 3 — S AT. 
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state any part of a quantum mechanical wave function 

that can be written in terms of exactly one of the basis 

elements \e n e n -i ■ ■ ■ ei). 
[3] In von Neumann's original approach the initial wave 

function is projected onto a full basis, where the choice of 

basis depends for the degenerate subspaces on the exact 

nature of the measurement apparatus. 
[4] Any choice of basis of the degenerate subspace leads to 

the same projection operator. 
[5] The application of Luders' postulate to entangled states 

is also of interest. 
[6] We initially rotate the first qubit to produce 000 • ■ • 01} + 

|000 • ■ ■ 00) and then step by step rotate the other qubits. 

This produces the desired wave function. 
[7] In the case of NP-complete problems, ala 3-SAT, instead 

of an oracle, one can associate the phase change to the 

solutions of the problem in hand. Each of the 2 N possible 

solutions has either its phase unchanged or gets a phase 

change of e tn . 

[8] Without loss of generality we set the marked state eu to 
be the fe-th record, i.e. e*., and equal to u + v, with u and 
v eigenstates of the first and second degenerate subspace 
of A respectively. Simple linear algebra shows that Cu = 
diu + d,2V with di = ai(l + it') and d-z = ^u'(ai + a 2 ), 



where u' is equal to etu and nonzero. Both d\ and di are 
nonzero. 

[9] If, for example, m is chosen to be log 2 (A r ) + 2 with 
log 2 N > 1, then the cumulative error probability is al- 
ways less then 1/3. 
[10] It could be that instead of the projection of the state onto 
a point in the degenerate subspace, one has a projection 
onto all the points in the degenerate subspace with an 
appropriate measure. This is not completely implausible 
to the author. 

[11] Besides the Wigner-Araki-Yanase theorem, what other 
restrictions are there on observables? How does one re- 
alize observables with degenerate spectrum? Are there 
construction mechanism, which are self-stabilizing? 

[12] An almost identical algorithm, to the one described in 
this paper, can find a marked record in a database with 
potentially a multiple number of marked records. Again 
the goal is to reduce, cycle by cycle, the subset of states 
to be considered, and eventually find at least one of the 
marked states. To be able to do this, one only has to 
remodel slightly Section llVl of the present algorithm. The 
essential fact that the degeneracy of the smaller subspace 
of A is broken in C remains unchanged. 



